Mon. Not. R. Astron. Soc. 000,[T]g](2013) Printed 20 January 2013 (MN MgX style file v2.2) 



Estimating the mass of the debris disc in HD 69830 

Kevin Heng 1 * 

1 Zwicky Fellow, ETH Zurich, Institute for Astronomy, Wolfgang-Pauli-Strasse 27, CH-8093, Zurich, Switzerland 



Submitted 2011 March 22. Re-submitted 2011 April 19. Accepted 2011 April 21. 



ABSTRACT 

We present a method to estimate the mass of the debris disc in the HD 69830 system, which 
also hosts three exoplanets with Neptune-like minimum masses. By considering the range of 
published stellar ages, we interpret the infrared emission from the debris disc as originating 
from a steady state, collisional cascade of dust grains. Using dynamical survival models sub- 
jected to observational constraints, we estimate the allowed range of disc masses. If the disc 
has an age rs 1 Gyr, then its mass is Afdi SC ~ 3-4 x 1O _3 M0, several times more massive 
than our asteroid belt. The maximum allowed age for the disc and the number of planetesimals 
it contains are determined by the assumed value for the binding energy of the planetesimals. 
If one insists on interpreting the disc as being transient, then this mass estimate becomes an 
upper limit. 
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1 INTRODUCTION 

Debris discs are believed to aris e from plane tesimals colliding and producing a cascade of dust grains, which reprocess the starlight into 



Debris discs are believed to arise Irom planetesimals colliding and producing a cascade ot dust grains, wnicn reprocess me starlignt into 
infrared emission dZuckerman 2001; Wvattl 2008 ). While infrared and sub-mm observations yield reasonable estimates for the dust mass, 



estimating the disc mass is generally challenging because the planetesimals cannot be detected via conventional methods jHeng & Tremaine 
2O10h . 

Of the hundreds of debris d iscs studied, the o ne residing in HD 69830 stands out because the system also hosts three exoplanets with 
Neptune-like minimum masses l Lovis et al. 2006h . Steady state models using a narrow interpretation of the disc properties infer that the 
infrared emission is anomalously bright, which implies that we may be observing HD 6983 at a special time in its evolu t ion — it is 
a "tran sient" debris disc fed by a recent, violent collision between planetesimals JWvatt et"af 2007 ; Heng & Tremaine 2O10h . Lisse et al. 
J2007ll reach a more subdued conclusion: based on the discovery of icy dust grains within the debris disc, they remark that the dust was either 
produced within the last year (relative to the time their observations were taken) or was being continuously replenished by planetesimal 
collisions. 

In this study, we offer an alternative, plausible perspective: as the stellar/disc age (assuming they are the same) is uncertain, it is possible 
to still interpret the disc emission as being non-transient. In this case, the old age of the disc becomes a friend and not a foe, because the 
possible values for the disc mass become severely restricted by both dynamical and observational constraints. In particular, we show that the 
disc in HD 69830 is several times more massive than the asteroid belt in our Solar System. In $2] we review the observational constraints on 
the HD 69830 system. In $3] we apply dynamical survival models to the system. We conclude by discussing the implications of our study in 
S|4] including the main model uncertainties. 



2 OBSERVATIONAL CONSTRAINTS 

The key aim is to meaningfully restrict the parameter space of our models using observational constraints. Therefore, we begin by 
reviewing the existing observational constraints found in the literature. Located 12.6 pc away, the HD 6 9830 star has a spectral type K0V, 
an effe ctive temperature T* — 5385 ± 20 K and a luminosity L* = 0.60 ± O.O3L0 dLovis et alJl2006h . By using theoretical evolutionary 
tracks. lLovis et al.1 $2W(k) derived a stellar mass of M, = 0.86 ± 0.03M Q . 



High precision radial velocity measurements using the HARPS spectrograph placed the exoplanets HD 69830b-d at semi-major axes 
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Figure 1. Left panel: dynamical survival models exploring the allowed range of disc masses as a function of the stellar age. Right panel: the corresponding 
allowed number of planetesimals. Both dynamically hot and collision-limited (denoted by "CL") discs are considered, with the latter being a special case of 
the former (see text). 



of a p — 0.0785, 0.186 and 0.630, respectively ( Lovis et al]|2006t) . The corresponding minimum masses are M p sini = 10.2, 11.8 and 18.1 
Earth masses, where i denotes the inclination of the system (assuming co-planar orbits). 

Infrared excesses at 12, 23.7 and 25 /im were detected in HD 69830; the optical depth (covering frac tion) of the associated dust grains 
was inferred tober = 3xl0~ 4 , more than 3 orders of magnitude higher than for the asteroid belt l Beichman et al. 20051) . The non- 



une 



detection of an infrared excess at 70 /im implies that the dust grains predominantly have sizes smal ler than 70 /xm/27r ~ 11 fj,m. The stron g 
mineralogical features detected in the spectrum from 8-35 /im further argues for small grain sizes (IBeichman et al 1120051 Jlisse et alj|2007h 
which rules out the possibility of the infrared excess being due to a dynamically warm disc (with particle sizes 2> 1 fim; Heng & Tremain 
2010). Detailed modelling of the infrared spectral energy distribution yielded semi-major ax es of a = 0.92-1.16 AU dLisse et al.11200' 
Constraints from mid-infrared interferometry place an upper limit of 2.4 AU on the disc size I Smith. Wvatt & Haniff 20091). These estimate s 
for a are consistent with stability analyses which suggest that planetesimals may exist at a = 0.8-1.2 AU l Lovis et al] 20061 : Ji et alj|2007 ). 

We thus adopt a conservative range of semi-major axes for the disc of a = 0.8-1.3 AU0 

The most uncertai n parameter in the s ystem is the stellar age t*. Ages of t* ~ 0.6-5 Gyr have been claimed (see Beichman et alliooj 
and references therein). IWvatt et al.l (12007) favour t* = 2 Gyr (see their Table 1); in tandem with assuming a single value of a — 1 AU, as 
well as fixed values for the parameters describing the plan etesimals, they conclude d that the infrared emission due to dust exceeds its steady 
state value, thereby qualifying it as a "transient" (see also lHeng & Tremainekoiol) . 



3 METHOD 



The formation of planetary systems from first principles is a challenging task fraught with theoretical uncertainties — these challenges 
are present for any formation model of the HD 69830 system ( Alibert et al. 20061) . Since our main interest is in estimating the disc mass, 
a simpler approach is to construct a survival model with the reasoning that any observed disc must have survived all dynamical processes, 
including gravitational instability, dy namical chaos, gravitati onal scattering, physical collisions and radiation forces that would lead to sig- 
nificant evolution over its lifetime I Heng & Tremaine 201oh . From a dynamical viewpoint, discs may be classified as being dynamically 
hot, warm or cold. Only dynamically hot discs are capable of producing a collisional cascade of dust grains. Collision-limited discs are a 
special sub-class of dynamically hot discs, where the collisional time scale is equal to the disc age. Dynamical survival models are relatively 
simple and may be described mainly by five parameters: the disc mass Afdisc, the disc semi-major axis a, planetesimal size r, planetesimal 
radial velocity dispersion a r and the s tellar/disc age t*. The ra dial velocity dispersion ay is in turn directly related to the root mean square 
eccentricity of the planetesimal orbits dHeng & Tremaing|20101) . 

Since observational constraints exist for a and i* (see ij2), we may thus use the models to constrain the allowed parameter space 
for the other three parameters. (The stellar mass M*, luminosity L+ and radius enter as minor parameters and are also fixed by the 
observations.) Additionally, an important constraint is provided by the observed dust optical depth r (see $2j- We permit only disc models 



1 Such a disc extends over one octave in semi-major axis, since it is centered at 1 AU and extends over Aa = 0.5 AU. 
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Figure 2. Masses of the exoplanets HD 69830b-d (dotted lines) as well as the spatial separation between them and the debris disc (in terms of Hill radii Rh; 
shaded regions), both as a function of the inclination angle i. 



with a theoretically calculated dust optical depth tq which satisfies the condition, 

e~ To < T < er . (1) 
The dimensionless quantity e allows for the uncertainties associated with our survival models, since the computed value of to depends on 
Mdisc r an d a r . Smaller values of e lead to narrower ranges of estimates for the disc mass. We choose to be conservative and adopt e = ylO 
such that the computed dust optical depth varies by a factor e 2 = 10. 

Figure Q] shows the range of possible disc masses Mdisc and the corresponding number of planetesimals in the disc TV as functions of 
t*. These diagrams were constructed by first considering a wide range of values for Mdisc, r and ay via a Monte Carlo approach. For each 
randomly generated trio of parameters, a list of dynamical conditions were checked, as described in detail in iHeng & Tremainel J201oh . The 
condition in equation l[T) was also checked. Only the discs which satisfy all of these conditions are retained. 

A few features of Figure Q~]are worth pointing out. Firstl y, Mdisc and TV are essentially unconstrained if the disc is allowed to be much 
younger than the star (e.g., as suggested bv lPavne et al .120091) . For young disc ages, the enormous allowed range of values for TV reflects the 
fact that a disc of mass Mdisc may be configured into a large variety of different planetesimal masses which still satisfy both the dynamical and 
observational constraints. In other words, the disc mass and observed infrared emission from dust grains are degenerate quantities. Secondly, 
if we allow the stellar/disc age to be t* ~ 1 Gyr, then decent estimates for the disc mass and number of planetesimals are produced: 
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Curiously, 

~ 10 _3 M 



Beichman et al 



(2005) extrapolated the inferred dust distribution to ~ 10-100 km-sized planetesimals and inferred a disc mass 
It is precisely because the disc is ol d and bright that th e range of possibilities for its mass is limited. 
It is important to note that the approach of Hen g & Tremaine J20ld> does not generally require TV 3> 1 for dynamically hot discs. 
For example, §3.2.1 of that paper describes a generalization of the criterion for the Toomre stability of a disc (which does assume a fluid 
treatment) to discs where TV ~ 1. In other words, while the dust grains are assumed to be in a collisional cascade, the planetesimals are not. 

It is reasonable to ask if any of the exoplanets in the HD 69830 system are capable of dynamically stirring the disc and thereby modifying 
the predicted dust co ntent away from its steady state value. Mutual gravitational interactions occ ur when two b odies are separated by less 



than about 3.5i? H Jchambers, Wetherill & BossI l996), where the mutual Hill radius is defined as d Ji et al ]|2007h 



Rh 



a + a p I M. 



3M* 



1/3 



(3) 



Figure |2] shows the spatial separations between HD 69830b-d and the disc, which we assume to extend from a — 0.8 AU to 1.3 AU. It is 
apparent that HD 69830b and HD 69830c are located at separations greatly exceeding 10T?h for non-zero values of the inclination. Only 
for inclinations of i < 4° is HD 69830d capable of dynamically stirring the disc, which corres ponds to Mp > 0.8Mj where M j is the 
mass of Jupiter. These values of i are not ruled out by stability analyses of the planetary orbits Lovis et al. 20061 : Ji et al. 2007 ). In such 
cases, the presence of HD 6983 0d may inc rease the radial velocity dispersion ay of the planetesimals in the disc and therefore increase the 
eccentricities of their orbits feafikovlbooil) . implying that c ollisions between them w ill become more violent. Non-intuitively, however, the 
collision rate either remains unchanged or decreases, since jHeng & Tremainell2O10h 

1k jconstant, G « 1, 

[a- 2 , e»i, 

depending on whether the self-gravity of the planetesimals is important (which is the case when the Safronov number <d exceeds unity). For 



© 2013 RAS, MNRAS 000,[T]g] 



4 Heng 



1-10 times the mass of Ceres, with radii ~ 1000 km, and 



the values of Mdi sc and N inferred from Figure Q] each of the planetesimals are < 
most likely have B > 1. 

Several mean motion resonances of the planetesimals with HD 69830d have already been identified bv lAlibert et al.l (120061) . Specifically, 
the 2: 1 resonance at a « 1 AU creates a zone of accumulation which allows the debris disc to survive. It is expected that the mean motion 
resonances will not substantially modify the dust production properties of the planetesimals away from the steady state model we have used 
in this study. 

We conclude that our dynamical survival models give reasonable estimates for the disc mass if i > 4° . 



4 DISCUSSION 



Using the existing observational constraints on the HD 69830 system, we have applied the dynamical survival models o flHeng & Tremaine 

(201(3) to estimate the mass of the debri s disc and find it to be several times more massive than the asteroid belt in our Solar System 
(~ 6 x 10~ 4 Mm; lKrasinskv et alj|2002l ). Our estimate hinges on being able to interpret the HD 69830 debris disc within the context of a 
steady state model, which is only possible if one adopts the lower range of stellar ages quoted for HD 69830 (i* « 1 Gyr). If one insists on 
inteipreting the disc as being transient, then this mass estimate becomes an upper limit. 

The most uncertain piece of physics involved in the dynamical survival models concerns the collisions between the parent planetesimals 
respo nsible for dust production. The binding energy per unit mass of a planetesimal is typically described by the function i Benz & Asphaus 
19991) . 



Q* = Qo 



m 
mo 



m 
mo 



(5) 



where the quantities Qo, mo, a and b are parameters determined from a combination of simulations and laborato ry experiments. The spec - 
ification o f the values for the indic es a and b also allows the mass distribution of the planetesimals to be fixed i Heng & Tremaine 2010). 
Following Heng & Tremaine 1 2010l) . we have adopted mo = 10 14 g, Qo = 6 x 10 5 erg g _1 , a = —0.13 and b = 0.44 by taking an average 
of the values assumed for ice and basalt, while being aware that they are some what uncertain. 

In particular, the value of Qo may vary by several orders of magnitude ( Stewart & Leinhardt 2009). For example, decreasing Qo to 
6 x 10 4 erg g _1 results in similar estimates for Mdi sc (equation |2j), but only allows for discs with t* < 0.7 Gyr to survive. The relative ease 
of breaking apart the planetesimals means that they are less likely to survive for a given time. If Qo = 6 x 10 3 erg g^ 1 , then we now have 
t* < 0.2 Gyr, which is inconsistent with the lowest value of the stellar age quoted in the literature (0.6 Gyr; see jQ . Also, as Qo decreases, 
the mass A/disc is distributed among a larger number N of planetesimals. For comparison, we note that lWvatt et al.l ( 120071) adopt a constant 
Q* = 2 x 10 6 erg g _1 . Our estimates highlight the fact that whether the disc is transient depends on the assumed value for the binding 
energy of the planetesimals. This quantity in turn depends on the material composition of the planetesimals, for which we only have Solar 
System asteroids and Kuiper Belt objects as a guide. 
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